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Abstract 

In the present paper we propose a doubly orbitally degenerate narrow-band 
model with correlated hopping. The peculiarity of the model is taking into 
account the matrix element of electron-electron interaction which describes 
intersite hoppings of electrons. In particular, this leads to the concentration 
dependence of the effective hopping integral. The cases of the strong and weak 
Hund's coupling are considered. By means of a generalized mean-field approx- 
imation the single-particle Green function and quasiparticle energy spectrum 
are calculated. Metal-insulator transition is studied in the model at different 
integer values of the electron concentration. With the help of the obtained 
energy spectrum we find energy gap width and criteria of metal-insulator 
transition. 
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I. INTRODUCTION 



The electron-electron interaction-driven metal-insulator transition (MIT) has fascinated 
theorists and experimentalists for many years. This transition is named after Sir Nevill Mott, 
being one of the pioneers who laid down the foundations of our physical understanding of 
this phenomenon [^0]. Despite a large amount of papers devoted to study of this transition 
(e.g. see a review [Q) the construction of consistent theory of MIT is still far from carrying 
out and constitutes one of the most challenging problems in condensed matter physics. 

The simplest model describing MIT in materials with narrow energy bands is the Hub- 
bard model 0]. This model describes a single non-degenerate band of electrons with the 
local Coulomb interaction. The model Hamiltonian contains two energy parameters: the 
matrix element to being the hopping integral of an electron from one site to another [to is 
not dependent on occupation of sites involved in the hopping process) and the parameter U 
being the intraatomic Coulomb repulsion of two electrons of the opposite spins. This model 
is studied intensively (for recent reviews see Refs. 0|1). 

Theoretical analyses, on the one hand, and available experimental data, on the other 
hand, point out the necessity of the Hubbard model generalization. As a rule this general- 
ization is perfomed by two means: taking into account the orbital degeneration or correlated 
hopping (in the present paper we do not consider the generalization of the Hubbard model 
by taking into account the inter-atomic Coulomb and exchange interactions). 

A model of non-degenerate band has to be generalized by taking into account the orbital 
degeneration being a characteristic of the transition metal compounds which exhibit MIT. 
First the degenerate Hubbard model was introduced for description of transition metal 
compounds in papers of Roth 0, Kugel' and Khomskii P], Cyrot and Lyon-Caen 0. In these 
works, in particular, the importance of intra-atomic exchange interaction J which stabilizes 
the localized magnetic moments in accordance with the Hund's rule was investigated. 

The intensive study of MIT in the degenerate Hubbard model has begun only in few 



recent years by use of the slave boson method [10-13|, the variational method WM, the limit 
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of infinite dimension . 



The results of works |T0|-p!5|] show the possibihty of the metal-insulator transition in the 
doubly degenerate Hubbard model. However the criteria of the transition obtained in the 
noted works are substantially different. In particular, at = 1, J = (here in all the 
cases T = K and the rectangular density of states are considered) in Refs. JTT| , P^ the 
criterion of MIT is U/w = 4.95, in Ref. ^ - U/w = 5.18, in Ref. ^ - U/w = 3.0; at 
n = 1, J = 0.1 in paper - U/w = 3.04; at n = 2, J = in Refs. [|ll|,|l2| - U/w = 6.0, 
in Ref. |10| U/w = 4.0, in Ref. lH - U/w = 9.0, in Ref. H - U/w = 3.7. In addition in 



paper ||T0[ at n = 1 MIT of first order is found, but in papers |jT2[ and [|Tj] - MIT of second 
order. Consequently a further investigation of MIT in narrow-band models with orbital 
degeneration is necessary. 

Another way of the Hubbard model generalization which allows to describe the essential 
peculiarities of transition metal compounds is taking into account the correlated hopping. 
The necessity of taking into consideration correlated hopping is caused by two reasons. 
Firstly, theoretical analysis points out the inapplicability of the Hubbard model for the 
description of real strongly correlated electron systems, in some compounds (e.g. see the 
estimation in Refs. |]I7|-pT|]) the matrix element of electron-electron interaction describing 
correlated hopping is the same order that the hopping integral or on-site Coulomb repulsion. 
Secondly, using the idea of correlated hopping and caused by it the electron-hole asymme- 
try we can interpret the peculiarities of some physical properties of narrow-band materials 

Now two ways are commonly used to generalize the single-band Hubbard model by taking 
into account correlated hopping. One of them has been proposed in Ref. [|16|. Hirsch showed 
that in contrast to the hopping integral of the Hubbard model (which is not dependent on 
occupation of sites involved in the hopping process) this integral of a generalized Hubbard 
model had to depend on occupation of sites involved in the hopping process. Hamiltonian 
of the generalized in such a way Hubbard model is written as 
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H = -Y. tlj^^ta^ja + UY^ rii^riii, (1.1) 

ija i 

tij = tAA(l - - + tAB{nia + rij^ - 2ni^nj^) + tBBniaUj^. (1.2) 



In recent few years Hamiltonian ( |1.1| ) is widely used to study MIT in narrow energy 
bands 

In Ref. P^p7[] the necessity of the Hubbard model generalization by taking into account 
the matrix element of electron-electron interaction describing intersite hoppings of electrons 
had been pointed out. The Hamiltonian of the generalized Hubbard model with correlated 
hopping is 

H = -fiY. (^ta^'i'^ + o.ta(^ja + {af^aj^Ui^ + /i.c.) + U rii^rin, (1.3) 

ia ija ija i 

with 

t{n) =to + nY,Jiikjk) (1.4) 

being the effective hopping integral of electrons between nearest-neighbor sites of lattice, 
X = J{inj)] n = Ne/N is the electron concentration [N^ is the number of electrons, is 
the number of lattice sites), 

2 

J^ikjk) = ip*{r - R,Mr - R^.)^-^|<^(r' - Rfc)|Mrdr', (1.5) 

ip{r — Rj) is the Wannier function, the prime at the sums signifies that i ^ j. 

In the model described by Hamiltonian ( |1.3| ) an electron hopping from one site to an- 
other is correlated both by the occupation of the sites involved in the hopping process (with 
the hopping integral X) and the occupation of the nearest-neighbor sites (with the matrix 
element J{ikjk) ai k ^ i, k ^ j) which we take into account in the Hartree-Fock approx- 
imation (Eq. (|1.4|) ). The peculiarity of model ( p..3|) is the concentration dependence of the 
hopping integral t{n) in contrast to similar models. 

MIT in a generalized Hubbard model with correlated hopping has been studied in a 



number of recent works p9|-p8|. In particular, at half-filling and to = ~X (or = 0) 



some exact results have been found [p9| , p2| -|3^ . In a simple cubic lattice with coordination 
number z MIT occurs at 



Uc = z{\tAA^tBB\) = 1z\tQ\. (1.6) 

If f/ > Uc the ground state of system is a paramagnetic Mott-Hubbard insulator with the 
concentration of doubly occupied sites (i = 0, the ground state energy is equal to zero. 

For an arbitrary to 7^ (or Iab 7^ 0) the finding of MIT criterion and the description 
of this phenomenon in a generalized Hubbard model with correlated hopping still remain an 
open problem. One of the step to solve this task is recent papers ||30| , |3ll , |35H 38|l where criteria 
of MIT, ground state energy, concentration of doubly occupied sites have been found. In 



Refs. ||3Qj3]]j35H 37| the authors have obtained the following criterion of MIT: 

Uc = z{\tAA\ + \tBB\) = z{\to\ + \to + 2X\) (1.7) 

in agreement with the Mott's general physical ideas 0. By means of the slave bosons 
method it has been found in Ref. that MIT occurs at Uc = ^z\t + X|; however here 
there is a problem of discrepancy of this result with the exact MIT criterion ( [1.61 ). 

Considering the above arguments on the necessity of the Hubbard model generalization 
by taking into account orbital degeneration, on the one hand, and correlated hopping, on 
the other hand, in the present paper we propose a doubly orbitally degenerate narrow-band 
model with correlated hopping. The structure of this paper is the following. Section H is 
devoted to the model formulation and model Hamiltonian construction, the representation 
of the Hamiltonian of a doubly orbitally degenerate model with correlated hopping in the 



electron and Hubbard operators is given. In Section |ITT| metal-insulator transition in the 
model at different integer values of the electron concentration is studied. The cases of the 
strong and weak Hund's coupling are considered. The absence of a natural small expansion 
parameter in the region near the MIT requires to find the nonperturbative approaches; in 
such a situation methods of mean-field type are useful. We shall use one of these methods. 



a variant of the generalized Hartree-Fock approximation |^,^ which has been proposed 



in Ref. [^,0. The approach gives the exact band and atomic hmits in the single-band 
Hubbard model and describes MIT. The method reproduces (see Ref. [^) the exact results 
(criterion of MIT (|1.6| ) and ground state energy) for the case of a half-filled non-degenerate 
band in a generalized Hubbard model with correlated hopping at to = —X. By means of the 
generalized mean-field approximation the single-particle Green function and quasiparticle 
energy spectrum are calculated. With the help of the obtained energy spectrum we find 



energy gap width and criteria of metal-insulator transition. Finally, Section |3 is devoted 



to the conclusions and the comparison of our results with those of other authors. 



II. MODEL HAMILTONIAN 



On the analogy of an orbitally non-degenerate model [|T3it3 start from the following 
generalization of the Hubbard model for an orbitally degenerate band taking into account 
the matrix element of electron-electron interaction which describes intersite hoppings of 
electrons (correlated hopping): 



+ U rii^^rii^i + U'Y^ Uiaani^a + {W - J)Y1 ^ia^^il^a, (2.1) 

where /i is the chemical potential, af^^, ai^a are the creation and destruction operators of an 
electron of spin a {a =1, |; a denotes spin projection which is opposite to a) on z-site and 
in orbital 7 (7 = denotes two possible values of orbital states), rii^a = o,tfa'^i-yo- is the 
number operator of electrons of spin a and in orbital 7 on i 

the hopping integral of an electron from 7-orbital of j-site to 7-orbital of i-site (we neglect 
the electron hoppings between a- and /3-orbitals). In real systems the electron hoppings 
between different orbitals can exist, in addition the hopping integrals are anisotropic for 
Cg orbitals. This may have an effect on the orbital and magnetic ordering [§,^,0]. We, 
however, assume for simplicity t'^j^ = tijSap- This assumption considerably simplifies the 
analysis of properties of the model under consideration and allows to describe the physics 
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e2 



of the metal-insulator transition. 

J{t^k^'nk^') = J J ip;{r - R,)cp^ir - R^.)^^:^|'^y (r' - Rfc)pdrdr' (2.2) 
{(p^ is the Wannier function), the prime at second sum in Eq. ( p.lD signifies that i ^ j, 

2 

U = l I |^,(r - R,)|2_l_|y,^(r' - R,)|2drdr' (2.3) 

is the intra-atomic Coulomb repulsion of two electrons of the opposite spins at the same 
orbital (we assume that it has the same value at a- and /3-orbitals), 

2 

U' = 1 1 - R^)\'J^^^\M^' - ROlMrdr' (2.4) 

is the intra-atomic Coulomb repulsion of two electrons of the opposite spins at the different 
orbitals, 

e2 



J = J J <(r - - R,)^-— ^v^;(r' - R,)<^„(r' - R,)drdr' (2.5) 

is the intra-atomic exchange interaction energy which stabilizes the Hund's states forming 



the atomic magnetic moments. The parameters U, U', J are connected by the relation |[45|| : 



U' = U - 2J. (2.6) 
In Hamiltonian (|2.1|) we rewrite the sum ' J(i'jk'j'j'jk'j')af^^nkYaj^^ in the form: 

(^J{i-fi'jj^i'j)a+^aj^^ni^^ + /i.e.) + (j{ililjlil)at^„aj^„ni^ + /i.e.) 
+ ' Y J{^lH 3lH)(^t^a(^h<^T^k^' (2.7) 

(7 = /3 if 7 = a, and 7 = a when 7 = /3). 

The first and third terms of Eq. ( p.7|) generalize correlated hopping introduced for 
an orbitally non-degenerate model (e.g., see Ref |2^). The second term of expression ( p.7|) 
describes correlated hopping of electrons being the pecuharity of orbitally degenerate systems 
only (it is absent in a single band case). Among this type of processes one can separate 



out three essentially different (from the energy point of view) hoppings (Xf '-representation 
allows to have done this easily). 

The first and second sums of Eq. ( p.7| ) describe the hoppings of electrons which are 
correlated by electron occupation of sites involved in the hopping process. The third sum 
describes the hoppings of an electron between |i70")- and |j7(T)-states which are dependent 
on the occupation number of other {k ^ i, k ^ j) sites. Let us take into account the 
influence of occupation of these sites in the Hartree-Fock approximation: 

J2 ' J2 J{hkj'jjk^')a+^aj^^nkY ~ nJ2'Ti{ij)a+^aj^„, (2.8) 

where n = {riia + ntp) is the average number of electrons per site, 

Ti{tj)=J2Ji^7k7nki) (2.9) 

(we suppose J^i'jkaj'yka) = J{i'ykPj'ykp) and Tiiij) have the same value for both a- and 
/?-orbitals). Assuming that a- and /5-states are equivalent, denote: 

J{tltljltl) = t'^M) = t'i3f,i^j) = t'ip (2.10) 
J{t^t^j^t^) = t'U^j) = t"pp{tj) = 4. (2.11) 

So we can rewrite Hamiltonian (|2.1|) in the following form: 

+ ^'{t'-jaty^aj.,„ni^^ + h.c.) + UJ2 '^nT^ni + U'Y1 ^iaanipa (2.12) 

ij'ycF i'y 



with the effective hopping integral tij{n) = tij + nTi{ij) being the concentration-dependent 
in consequence of taking into account correlated hopping Ti{ij). 

The distinction of Hamiltonian ( p^.l2D from models of narrow-band materials with orbital 



degeneracy is taking into account the matrix element J{i'yk'y'j'~fk'y') caused by electron- 
electron interaction. This leads to the electron-hole asymmetry (which is analogous one to 
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the case of an non-degenerate band [^,^) and the dependence of hopping integral on the 



average number of electrons per site. Thus the narrow-band model which are described by 



Hamiltonian (|2.12| ) shows much better physics than the Hubbard model with doubly orbital 
degeneration. 

In the model described by Hamiltonian ( |2.12| ) each site of the lattice can be in one of 16 
states (see Fig. p. 



Let us pass to the configurational representation of Hamiltonian (|2.12|) . Use the repre- 
sentation of the operators of creation and destruction of electron through '-operators of 
transition of site i from the state / to the state k: 



(2.13) 









+ Xj^'^^ 






+ Xf Tar 




+ ^4,,2i 








+ X/T-'^T 


-tAa2.Q!t 












= xr 


_ ji^TT.oT 


— X^^''^^ 


+ Xf ^'^^ 






+ Xf^^'"^ 


+xr^ 
















+ Xf ^'^T 


_ j^4,a2T 



Such choice of representation with taking into account the commutational rules 

{Xf''; Xj'^} = Sij{SikXf'^ + StpXf'^), 

[Xf; Xf ] = 5,:,(5,,Xf'* - 5,,Xf (2.14) 

where the anticommutator {{A; B}) has to be taken only if both operators are fermionic, 
i.e., change the particle number by one (e.g. X7'^'° or X/^'^'^), and the constraint 

E^f = l, (2-15) 

p 

where Xf = Xf''xj'^ is the operator of number of \p) states on the site i ensures the 
fulfillment of the anticommutation relations for a-operators. 

The model Hamiltonian in the configurational representation has the form: 

+(f/' - J) E xr + f/' E + E xf (2-16) 



la la %^( 

-{U + 2U' -J)Y. xf" + 2{U + 2U' -J)Y. Xf + Ht, 

i'ya i 



where the kinetic part of the Hamiltonian is 

Ht = ^ Hnm (2-17) 



n,m 



with n, m = {0-7(7, -ya-aa, -ya-aa, 7(7-72, aa-j2a, aa-'~^2a, ^2-^2a, 72(7-4} . 

The Hamiltonians Hnm contain the processes which form the energy subbands (analogues 
of the Hubbard subbands) and the processes of the hybridization of these subbands. In Fig. |^ 
the transitions between the states of the sites which form the corresponding subbands are 
shown. 

The different hopping integrals t^™ correspond to the transitions within the different 
subbands if„„ or between the different subbands Hnm {n ^ m): 

nr = U, + {T^ + T"'), (2.18) 

where 



^0-70- _ 


= 0, 




-72(T 


-t'--+ t"- 


^70--CT(T 






-■y2a 


-t'-- + t"- 




- 1' 


^72- 


■720- 


= Kp 


^7o--72 


- 1" 


7-72<T 


-4 _ 





(2.19) 



The mutual placement and the overlapping of the noted subbands depend on the rela- 
tions between the values of intra-atomic Coulomb repulsion parameters U, U', intra-atomic 
exchange interaction parameter J and the widths of subbands. 

III. METAL-INSULATOR TRANSITION 

At integer values of electron concentration n = 1, 2, 3 in the model described by the 
Hamiltonian ( |2.16| ) MIT can occur. The possible metal-insulator transitions at some integer 
values of the mean electron number per site will be considered below. 
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A. Metal-insulator transition at electron concentration n = I 



1. The strong Hund's coupling case 

Let us consider the case of the strong intra-atomic Coulomb interaction U' ^ tij and 
the strong Hund's coupling U' ^ U' — J (values U' and J are of the same order). These 
conditions allow us to neglect the states of site when there are more than two electrons on 
the site and the "non-Hund's" doubly occupied states | ti), | it), 1^2), |/92) (the analogous 
conditions are used for an investigation of magnetic properties of the Hubbard model with 
twofold orbital degeneration in Ref. [^5|-47|). Thus lattice sites can be in one of seven 



possible states: a hole (a non-occupied by electron site) |0); a single occupied by electron 
site lat), l/^T), the Hund's doublon or site with two electrons on the different 

orbitals with the same spins | ft)) I ii)- 

Let us pass to the configurational representation of the Hamiltonian. Electron operators 
in terms of the transition operators Xf'' of site i from the state |/) to the state \k) are 
expressed by the formulae: 



(3.1) 







+xr^^\ 






+ Xf 


+ 








_ -^0,01 






= xr 






= xr 






= xr 




«j/3i 


= xr 





(in these formulae we have took into account only seven possible states). It should be 
emphasized that one have to consider all the possible site states in order to get the correct 
anticommutation rules for electron operators. 

Let us calculate quasiparticle energy spectrum of the model neglecting the correlated 
hopping t'ij = t'-j = 0, tij{n) = tij. We write the Hamiltonian in the form: 

H = -/i^ (xr + Xf" + 2Xr) + (X7"'°X°'^" + Xf^'^'X/"'"") + 

icr 

+ (U'-J)J2xr + H^p, (3.2) 
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where 



X^^'^^Xf^^ + h.c.). (3.3) 

For calculation of quasiparticle energy spectrum and polar state concentration c = (Xf) 
- holes, = {XP), di = (X/^) being the Hund's doublons we use the method of the 
retarded Green functions. For the Green function 



(3.4) 



we have the equation: 



(i?-/i)((x,"^X'"')) 



(3.5) 



where 



(3.6) 



tjja 



To break off the sequence of equations for Green functions we apply a variant |42] of the 



generalized mean field approximation |40,41 



Xf^'^,Ho\=Y.'<Pj)X- 



Xf'',H^^]=J2'e,ipj)X] 



TT,/3T 



(3.7) 



where €{pj), ei(pj) are the non-operator expressions, the method of its calculation will be 
given below. The representation ( p.7| ) of commutators in Eq. ( p.5[ ) is defined by the operator 
structure of Hq and Hais. 

Using the approximation ( |3.7| ) Eq. ( |3.5| ) for the Green function G^pp,{E) can be written 

as 



Srm' 



(E - n){{xfV/^)) = ^(x; + xf) + T.'<pj){{xf'V/')) 



(3.8) 
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For the Green function 



we write the analogous equation: 

3 

+E'^2(pj)((x;T''^;;"^)), 

3 

where the non-operator expressions e, €2 are defined by the formulae: 



(3.9) 



(3.10) 



3 



(3.11) 



Thus, we obtain the closed system of equations for functions G'^^1{E) and G^^ME). 



.(2), 



To calculate e(pj), ei(pj), e(pj), e2{pj) we use the procedure proposed in paper p2 
The values of e{pj), ei{pj), €{pj), €2{pj) we find by anticommutation of Eqs 



with the operators X^;"^ and X^,^'^^ respectively: 



(Xj + Xf)e{pp') = {x'/\ [Xf'', H,] } , 
(X^ + X^)e,{pp') = {xf;'^T. [x;T,o^ H^^] } , 
+ Xl^ipp') = {X^^/'^T. [^n,m, } , 



(3.12) 



We rewrite Xf '-operator in the form Xf'' = a^^an, where a^, an are the operators of 
creation and destruction for \k)- and |Z)-states on z-site respectively (with the constraint 



Y.c4kCiik = 1); thus Xf = ataio, X 



(T- (T- X 



7fT 



aj^^aj^o-- Let us substitute a- 



operators by c- numbers (here there is a partial equivalence with the slave boson method [^ ). 
At n = 1 in a paramagnetic state 



d = d^ = d^ = c/2, (Xf ^) = 0.25 - d 



a. 



1/2. 



a, 



Oil 



a. 



i'ycr 



ai-y„ = (0.25 — d) 



1/2 



(3.13) 
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The proposed approximation is based on the following physical idea. When the lower 
and upper Hubbard subbands overlap weakly (the case of a paramagnetic Mott-Hubbard 
semimetal) |0)- and |crcr)- states almost do not influence on |7cr)-states then 1 70") -subsystem 
can be considered as a quasiclassical subsystem (an analogue of the thermodynamic reser- 
voir). In this case one can substitute a-operators through c- numbers. We extend this 
ideology to determination of e{pj), ei(pj), e{pj), e2{pj)- Thus we substitute the creation 
and destruction operators of |0)-, \<J<j)- and |7cr)-states through respective quasiclassical ex- 
pressions. Actually the proposed procedure is equivalent to a separation of the charge and 
spin degrees of freedom. 

Then we obtain the formulae in k-representation: 

„ , -O.Sc^ + 0.5c -0.25 „ , 2n 

= 0,25 + 0.5c = ^l'"' + 

E(k) = (c^ + c - 0.5)i„ ,M = °-^^^^t^. (3.14) 



From system of Eqs. ( |3.8|) , (|3.10|) we have 



(1) 0.5c + 0.25 E-^^ + U'-J-~eik) 

~ ^ {E + EM){E + E,{k)y ^^-^^^ 

r^(2).p^_ 0-5c + 0.25 .,(k) 

- 2n {E + EM)iE + E,ik)y ^^■^^> 

where the expressions 

EM - - + i{Kk) + .(k) - (£/' - Jff- 

+4e(k)((7' -J)+ 4l£i(k)e2(k) - e(k)e(k)|}"' (3.17) 

are the quasiparticle energy spectrum. E2{k) describes the electron spectrum in the 0-aa- 
subband (an analogue of the lower Hubbard subband), Ei(k) describes the electron spectrum 
in the aa-crcr-subband (an analogue of the upper Hubbard subband) . The computation of the 
quasiparticle energy spectrum in the O-^a- and /Ja-acr-subbands gives the same Eqs. ( ^.171 ). 

The energy gap width (difference of energies between bottom of the upper and top of 
the lower bands) is 
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AE = E2{-w) - Ei{w). (3.18) 

The peculiarity of the obtained energy spectrum and energy gap width is their depen- 
dence on the hole concentration c. With rise of temperature or the parameter [U — J)/ (2w) 
the holes concentration decrease smoothly. 

The energy gap width AE as a function of the parameters {U' — J)/ (2w) and (kT) / {2w) 
is presented in Fig. § and Fig. ^ respectively. With a change of the parameter {[/' — J)/ (2w) 
system undergoes the transition from an insulating to a metallic state (negative values of 
the energy gap width correspond to the overlapping of the Hubbard subbands). In the 
model under consideration at T = K insulator-metal transition at n = 1 occurs when 
{U' — J)/{2w) = 0.75 (Fig. 1^, the lower curve) (in the single-band Hubbard model the 



respective parameter is f// {2w) = 1 p2 



The transition from a metallic to an insulating state with increase of temperature at 
given value of the parameter {W — J)/{2w) is also possible (Fig. It can be explained by 
the fact that energy gap width AE ( |3.18| ) increases at increasing temperature T which is 
caused by the rise of polar states concentration at constant w, {If — J). 

2. The limit of the weak Hund's coupling 

Let us consider the MIT at electron concentration n = 1. If the exchange interaction 
is small comparatively to the Coulomb interaction J <^ U then we can take J into ac- 
count in the mean- field approximation (see, e.g., Ref. Ijl0|)- The corresponding terms of the 
Hamiltonian ( p.l2|) can be written as: 



-2 J ^ rtiaariipa " 3 J ^ UiaaUifSa = "2 J '^{{niaa)nii3a 
ia icr ia 

+niaa{nii3a)) " 3J '^{{niaa)nif3a + niaa{nii3a)) = {3.19) 
ia 

5t/ ("^ij-yo-) ^ ^ ^i-ya- 



i^a 

The Hamiltonian ( |2.12[ ) takes the form: 
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i'ycr ii70" 

~l~ ^ ^ (tij^i'^/cr^j'yo-^i^ ~^ h.C.) + ^ ^ {t^jd^^^Clj^^jTli^^ + h.C.) (3.20) 
~\~U ('^27cr^27f7 ~t~ '^i'ya^i^a ~t~ '^i'ycr'^i'yo'^ ^ 



icT'y 

where ft = fi + 5J(nj^o-). 

Considering MIT at the electron concentration n we can take into account in the Hamil- 
tonian only the states of site with n — 1, n, n + 1 electrons (the analogous simplification 



has been used in Refs. |[TT|JT3|1 ). In the vicinity of the transition point at the electron con- 
centration n = 1 the concentrations of sites occupied by three and four electrons are small. 
Neglecting the small amounts of these sites we can write the electron operators in the form: 



(3.21) 





_ j^a],0 




+ xj^^f"^ 


1 va2,al 
+ ; 






_|_ J^ii:/3i 


+ XP'^^ 










— X^'^'"'^ 


+ Xf '^^ 


^tf3i 


= xr 


— x^^'"^ 







Hamiltonian ( |3.20D in the configurational representation takes the form: 



+ [/ I 5: Xr + E + E + H2, (3.22) 



7 la t-y 

where the kinetic part of the Hamiltonian is 



^1 = E (t^^ E xr'xp^ + (t,, + 2t'r) E xr^'^x 

i,j \ 7°" 7°" 

+{tij + 2ty (^^xr^'x;'^'"" + E(^4"'''""^r'''' 

7cr 



7(7,72 

i 



+Xf '^^Xf + (Xr'""X""'''" + Xr^'Xf + /i.e.)) ) (3.23) 



70" 

+ E vAxr^'^'xp''^ + xp'^^xp'"^) + h.c. 

cr 
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ij \ -yo- 

- Xf"'°X;^'^" + h.c.)) (3.24) 

+{u,+t':,)T.iv.xr'xp^'+h.c.)], 

ja J 

where rji = rja = 1, Vi = = 

The processes that form energy subbands are included in the Hamiltonian Hi, the pro- 
cesses of hybridization of these subbands are included in the Hamiltonian H2. 

The single-particle Green function can be written as: 



(3.25) 



Here we have introduced the following notation: Yp = X^'°^ + X^'^ +X^^''^K The functions 
((X0'"TpCp1'°)) and ((XO'^TfKpt)) satisfy the equations: 



(X° + XpJ) 



{E + f,- U){{Yp\X;}'')) = {{[Yp, Hi + H,]\X^'^')). 



(3.26) 



For calculation of these functions we use the generalized mean-field approximation Let 
us take the commutators in (|3.26| ) in the form: 



[Yp, Hi] = E'C(PJ)^., [Y,, m = J2'aPj)X^'"\ (3.27) 
j j 

where e{pj),e{pj),({pj),({pj) are the non-operator expressions. 

After transition to k-representation the system of equations ( |3.2(j| ) taking into account 

( 3.271 ) has the solutions: 



{{Y,\X^}'')h 



27r {E - Ei{k)){E - E2{k))' 
(X° + Xf) C(k) 

27r {E - Ei(k)){E - E2{k))' 
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(3.28) 



Here £'i_2(k) is the quasiparticle energy spectrum: 



, , U , e(k) + C(k) ^l 



T ^ {[f/ - e(k) + C(k)]2 + 46(k)C(k)}'/' . (3.29) 



The same expression for the quasiparticle energy spectrum we obtain from the functions 
((Xp'"1^|yp"))k and ((y^|ypt))k- The values of e{pj),i{pj)X{pj), C{pj) "we find by anticommu- 



tation of (|3.27| ) with the operators X"J' and Ft: 



{X^ + X^ + + XtT + Xj + Xfre{pp') = {y;- [x^°'-T^ } ^ (3.30) 

(x^ + + + x^; + xji + xfxipp') = {y;t; [f„ /Ji]} , 

By use of the mean-field approximation analogously to the above, in the case of = t'l we 



obtain 



e(k) = (c + fe) + 



3b^ ' 
c + b. 



3tk Ti 

c + 6' 



d + h d + bj 



(3.31) 



C(k) = f + 



52 6d'' 

+ 



2 c + b c + b J 

here c, 6, d are the concentrations of the holes and sites occupied by one, two electrons, 
respectively, connected by the relations: 



c = 6d, b = - — 3d] 



(3.32) 



and 



tk — t]i + 2t']^. 



(3.33) 



In the point of transition, when the concentrations of the holes and doublons are equal to 
zero, the energies of the electrons within the subbands are 



E2{k) = -jl + U + (3.34) 

The energy gap in this case is 

AE = U-w-w = 0, (3.35) 
where w = z\tij\, w = z\tij\. From the equation ( p. 35 ) we obtain the criterion of MIT: 



U = w + w. (3.36) 
In the partial case = = (i^ ^^is case = i^) we have 

In the Fig. |^ the energy gap for different values of the correlated hopping at T = is plotted. 
With the increase of the correlated hopping at the fixed value of parameter U / 2w the energy 
gap width increases and the region of values of U/2w at which the system is in a metallic 
state, decreases. 

B. The limit of the weak Hund's coupling at electron concentration n = 2 

Let us consider the MIT at electron concentration n = 2. In the vicinity of the transition 
point in the case of two electrons per atom the concentrations of holes and sites occupied 
by four electrons are small. Neglecting the small amounts of these sites we can write the 
electron operators in the form: 

n+ - vUM , vP'iM x-/32i,ii , x-a2T,a2 

^+ _ Yn,al Yl32,(i1 , Ya2l,a2 , v./32T,TT , yI321,U 
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Let us rewrite the Hamiltonian ( |3.2CI| ) in the configurational representation at electron con- 



centration n = 2. For the small values of the intra-atomic exchange interaction (J ^ U) 
we take J into account in the mean- field approximation ( |3.19| ). The Hamiltonian takes the 
form: 



\ 270" i(j 27 270" 

+ u\y^ xr + E + E ^7' + 3 E +h, + h„ (3.39) 

\ ia ia i'y i'ya j 

where the kinetic part of the Hamiltonian is 



^ 70" cr 

+ + xP'^^Xf"'"" + h.c.)) ] (3.40) 

+(i..+4+i:-)(E^<^^7^r"^^"'' 

7cr 

+ + 2t;;.) (E x7''^'"'^xf '^'-^ + E(^^°'"'""^i 

70" cr 



7cr 

+ + 2t^ + 2t;'.) ^(xf -.-xf '"'-^ + xf -'--xf ''^'^ + /i.e.) 



7cr 

+ 3t^ + t%)Y^^,{Xf''~^''xf^^'"' - xf '^''^^Xf '"''^ + h.c) 



h2 = y: (iu, + 2t^ + 1:^) (Y^ivaxr^'-xp'''' + n-^xf^^^^xf^^^-^] 

ij \ ^7^ 
a 

H-Aj Aj- j+r^^jl^Aj A^- -Aj Aj- j 

+(^.. + 3t^)(E^^^r'^''^i''^'" + Y.ixf'^'^xp^^'' 
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70" 

The processes that form energy subbands are included in the Hamiltonian Hi, the processes 
of hybridization of these subbands are included in the Hamiltonian H2- 
Let us write the single-particle Green function as 

(W«tI4.t)) = ((^P^P' )) + ((^pF^pt)) + mz;,)) + {{Z,\Z;,)). (3.41) 

Here the following notations have been introduced: Yp = X^'"^ + X^'^ + X^^'^^, Zp = 
Xp^'"2i + XpiT'"2T +x^2,/32T^ rj^Yie functions {{Yp^pt)) and {{Zp^pt)) satisfy the equations: 

{E + fi- U){{Yprpt)) = ^S,j + {{[Yp; Hi + H,Wp^)), 

(E + f^- 2f/)((Z,^+)) = {{[Zp; Hi + H^Wpt)), (3.42) 

where A = {X^ + X^ + X^^ + X]} + Xji + Xf). On the analogy of the previous section 
we use the generalized mean-field approximation to calculate these functions. Let us take 
the commutators in ( p.42| ) in the form: 



[Yp, Hi] = Y!<P3)Yv [Yp^ H,] = Y!~<V3)Z,, 

j j 

[Zp, Hi] = E'C(PJ)^., [Zp, H,] = Y.'CiPj)Yv (3-43) 
j j 

where e{pj),e{pj)X{pj)/C{pj) are the non-operator expressions, which we calculate using 
method of paper |P2| . 

After transition to k-representation the system of equations ( |3.42| ) taking into account 
( p.43| ) has the solutions: 



((YY^)) E + f^-2U-C{k) 

U^pIV //k (E - Ei{k)){E - E2{k)) ' 

^^'^^^'^^^ - T^ iE-EiimE-EMY ^'-''^ 



Here ii^i^2(k) is the quasiparticle energy spectrum: 
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= + f + ^^^^^ t\{[U- 6(k) + C(k)]^ + 4e(k)C(k)}^/^ . (3.45) 

The same expression for the quasiparticle energy spectrum we obtain from the functions 
{{Yp\^p,))\^ and ((Zp|Z^))k. The values of e{pj),e{pj),({pj)X{pj) we find by anticommuta- 
tion of (|3.43|) with the operators l^t and Zp,: 



{X^ + + + X^J + + Xf)e{pp') = {y;-, [Y„ H,] 
{Xj^ + Xj;^ + Xf + Xf T + xf i + Xf T)g(pp') = [Y„ H,]} , 

(X^ii + + xf + Xf T + xf i + Xf T)C(pp') = [Z„ /Ji]} , (3.46) 

By use of the mean-field approximation analogously to the above, in the case of = t'^ we 
obtain 

e(k) = tk 3(rf + b) + - — - + — — - - tl- 



d + b 3{d + b)J "Sid + b)' 

,-(k) = ^f3(d + 6)-^ + ^V (3.47) 

^ ' 2 V ' 3{d + b) 3{d + b)J' ^ ' 

C(k) = -4:;7l^ + tkf3(d + 6) + -^+ ^ 



^3{d + b) V ' d + b 3{d + b) 



with 



3(d + 6) 3{d + b)J' 
= t]i + 2t'^. 

tt = tk + 44; (3.48) 

here b is the concentration of the sites occupied by one (or three) electrons, d is the concen- 
tration of the doubly occupied sites, connected by the relation: 

l-8d , , 

b = 3.49 

In the point of transition, when the concentrations of the singly and triply occupied sites 
are equal to zero, the quasiparticle energy spectrum is 
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Using the quasiparticle energy spectrum ( p.50| ) we find the energy gap width. In the 
point of MIT the energy gap is equal to zero. From this condition we find the criterion of 
MIT. In the partial case of t'^ = = (in this case tj. = tk) we find 

^ = (3.51) 

In the Fig. |^ the energy gap for different values of the correlated hopping at T = is 
plotted. With the increase of the correlated hopping at the fixed value of parameter U/2w 
the energy gap width increases faster than at = 1 and the region of values of U /2w at 
which the system is in the metallic state, decreases, analogously to the case n = 1. 

C. The limit of the weak Hund's coupling at electron concentration n = 3 

Let us consider the MIT at electron concentration n = 3. In the vicinity of the transition 
point in the case of three electrons per atom the concentrations of holes and sites occupied 
by one electron are small. Neglecting the small amounts of these sites we can write the 
electron operators in the form: 



(3.52) 





















_ J^4,/32T 












« = 


_|_^a2i,a2 


_^ ^/32T,TT 




_ j5^4,a2T 



Let us rewrite the Hamiltonian (|3.2CI| ) in the configurational representation at electron 



concentration n = 3. For the small values of the intra-atomic exchange interaction (J ^ U) 
we take J into account in the mean field approximation ( p. 19] ). The Hamiltonian takes the 
form: 
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H=-ii[2j2{xr+xn+2j: xr + 3 xr^^ + 4 ^ 

\ 1(7 i-y i'ya i 

+ uIy: Xr + E + E ^7' + 3 E Xf" + &Xt\+H, + H2, (3.53) 



where the kinetic part of the Hamiltonian 



a 

= Yi^k, + 3t^ + 24) E(^7^7"^"'^^r'' 

+(^.. + 3t:,) E^<^^7'""'^r'' + ■ 

70- / 
The processes that form energy subbands are included in the Hamiltonian ifi, the processes 
of hybridization of these subbands are included in the Hamiltonian H2- 
Let us write the single-particle Green function as: 

((«paTl4aT)) = ^kXf^'^P'')) + ((^p'^'^p^)) (3-54) 

Here the following notations have been introduced: = X^'^^i +X^T,«2T ^ j^^2,^2T ^^-^^ 
functions ((X^2i4p^4,M^^ l^^^^i^f^ satisfy the equations: 

+ ^ - 3f/)((x,^^^p^;;^^)) = + mf"': Hi + H,]\x'/')), 

iE + il-2U){{Z,\X'/'')) = {{[Z,;H^ + H,]\X^^'')). (3.55) 
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On the analogy of the previous sections we use the generahzed mean-field approximation to 
calculate these functions. Let us take the commutators in (|3.55| ) in the form: 



j j 



(3.56) 



where e{pj),^{pj),({pj)j C{pj) are the non-operator expressions. 

After transition to k-representation the system of equations ( |3.55| ) taking into account 
( p.56| ) has the solutions: 

(X^'^ + X^) E + fi-2U-C{k) 



27r {E - E^{k)){E - E^ik))' 
(^p^^ + ^p^) C(k) 



{{ZplX^'^))^ 27r (E - E,{k)){E - E2{k))- 

Here ii^i 2(k) is the quasiparticle energy spectrum: 

i^i,2(k) = -A + ^ + ^i^l±^ tU[U- e(k) + C(k)]^ + 4e(k)C(k)} 



(3.57) 



1/2 



(3.58) 



The same expression for the quasiparticle energy spectrum we obtain from the functions 

The values of ^{pj),e{pj),({pj), C{pj) we find by anticommutation of ( |3.56| ) with the 
operators X^;^"^ and Zi: 



(3.59) 



ix^'' + K'Hpp') 

{X^ + X^ + Xf + Xf T + xf i + Xf T)e(pp') 
{x^ + x^; + Xf + Xf T + xf^ + Xf t)c(V) 

{X^,'' + Xl)~app') 



By use of the mean-field approximation analogously to the above, in the case t'^ = t'l we 
obtain 

^(k) =t'Jit + f) + ^]- 3t* 



{X'/'';[X^''\H,]} 
{Z+;[X;'-T,ij2]}, 
{Zp'] [Zp,Hi]}, 
|Xp,^'°; [Zp, H2 
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,„.!:^(,,;,_JL,_^), (,ao) 

C(k) = -tl-L- + tl 3{d + t) + -—+ ^ 



'■d + t ' d + t (d + t) 



here 



t + f t + fj' 
= tk + 4t^, 

tk = + 64; (3.61) 

d, t, f are the concentrations of the sites occupied by two, three and four electrons, respec- 
tively, connected by the relations: 

f = 6d, t=^-3d. (3.62) 

In the point of transition, when the concentrations of the holes and single electrons are 
equal to zero, the energies of the electrons within the subbands are 

Ei(k) = -jl + 2U + tl 

E2{k) = -fl + 3U + tl. (3.63) 

The energy gap in this case is 

AE = U - w* - w' = 0, (3.64) 

where w* = z\t*A, w* = z\t'A. 



From the equation ( |3.64|) we obtain the criterion of the MIT at the electron concentration 
n = 3: 

U=w* + w\ (3.65) 
In the partial case 4 = 4 = (in this case tk = tk) we have 

— = 1. (3.66) 
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This result coincides with the corresponding critical value at the electron concentration 
= 1 in the consequence of the electron-hole sjTiimetry of the model without the correlated 
hopping. 

In Fig. 1^ the energy gap for different values of the correlated hopping at T = is plotted. 
With the increase of the correlated hopping at the fixed value of parameter U / 2w the energy 
gap width increases faster than at n = l,n = 2 and the region of values of U/2w at which 
the system is in a metallic state, decreases. 

IV. DISCUSSIONS AND CONCLUSIONS 

In the present paper we have proposed a doubly orbitally degenerate narrow-band model 
with correlated hopping. The peculiarity of the model is taking into account the matrix 
element of electron-electron interaction which describes intersite hoppings of electrons. In 
particular, this leads to the concentration dependence of the hopping integrals. Using the 
representation of Hamiltonian of a doubly orbitally degenerate model with correlated hop- 
ping in terms of the Hubbard operators the cases of the strong and weak Hund's coupling 
have been considered. By means of a generalized mean-field approximation we have calcu- 
lated the single-particle Green function and quasiparticle energy spectrum. Metal-insulator 
transition has been studied in the model at different integer values of the electron concen- 
tration. With the help of the obtained energy spectrum we have found energy gap width 
and criteria of metal-insulator transition. 

The peculiarities of the expressions for quasiparticle energy spectrum and energy gap are 
dependences on the concentration of polar states (holes, doublons at n = 1; single electron 
and triple occupied sites at n = 2; doublons and sites occupied by four electrons at ri = 3), on 
the hopping integrals (thus on external pressure). At given values of U and hopping integrals 
(constant external pressure) the concentration dependence of AE allows to study MIT under 
the action of external influences. In particular, Ai?(T)-dependence can lead to the transition 
from a metallic state to an insulating state with the increase of temperature (see Fig. the 
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described transition is observed, in particular, in the (Vi_KCr^)203 compound 0,0 and 



the NiS2-a;Se2^ system P9| , |50| . The similar dependence of energy gap width can be observed 
at change of the polar states concentration under the action of photoeffect or magnetic field. 
The strong magnetic field can lead, for example, to the decrease of polar state concentration 
(see Ref. |^B|) initiating the transition from a paramagnetic insulator state to a paramagnetic 
metal state. Contrariwise, the increase of polar state concentration under the action of 
light stimulates the metal-insulator transition, analogously to the influence of temperature 
change. At the increase of bandwidth (for example, under the action of external pressure or 
composition changes) the insulator-to-metal transition can occur. 

The results allow to study the influence of the correlated hopping and orbital degeneracy 
on MIT. The dependences of energy gap width on the parameter U /2w in absence of the 
correlated hopping (tj^ = = 0) different electron concentrations are given in Fig. p. 
One can see that in the case n = 2 the MIT occurs at smaller value of U/2w then at 
n = 1. This result is in qualitative accordance with the results of work [ffO], in distinction 



from [^,^. Using the critical values of the papameter U/{2w) at which MIT occurs for 
different integer electron concentrations (see Fig. |^) we can interpret the fact that in the 
series of disulphides MS2 the C0S2 (one electron within Cg band corresponding to n = 1) and 
CuS2 compounds (three electrons within e^-band corresponding n = 3) are metals, and the 
NiS2 compound (two electrons within e^-band corresponding n = 2) is an insulator. Really, 
for 0.94 < U/2w < 1 at the electron concentration n = 2 system described by the present 
model is an insulator, whereas for the same values of the parameter U/2w at the electron 



concentrations n = 1, 3 system is a metal (according with the calculations of Ref. the 
ratios U/2w in these compounds have close values) . 

We have found that in the case of the strong Hund's coupling at n = 1 metal-insulator 
transition occurs at smaller value of the parameter {{U — J)/2w)c = 0.75 than in the case 
of the weak Hund's coupling (([/ — J)/2w)c = 1. 

When the magnetically ordered states are taken into account the phase diagram of the 
considered model (Fig. |^) can be changed. In particular, with increase of correlation strength 
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the transition from paramagnetic to magnetically ordered state (antiferromagnetic in- 
sulator or ferromagnetic insulator) can occur, similarly to the magnetic transition found in 
the Ref. |]T2[ by use of slave-boson method for the doubly degenerate Hubbard model. 

At nonzero values of correlated hopping the point of MIT moves towards the values of 
parameter U /2w (Figs. |^-|^) at which system is a metal in proportion to correlated hopping 
value (Fig. |10D. From Fig. |lOl one can see that {U/2w)c decreases with increasing correlated 
hopping, and what is more at n = 2 with the increase of the correlated hopping parameter 
K the value {U/2w)c decreases faster than at the electron concentration n = 1, and at n = 3 
the value {U/2w)c decreases faster than at n = 1, n = 2. The non-equivalence of the cases 
n = 1 and n = 3 is a manifestation of the electron-hole asymmetry which is a characteristic 
of the models with correlated hopping. 

Thus both orbital degeneracy and correlated hopping are the factors favoring the transi- 
tion of system to an insulating state in the case of half-filling with the increase of intra-atomic 
Coulomb repulsion in comparison with the single-band Hubbard model (in this connection 
see Refs. ^J^). 

In the present paper considering MIT we have neglected the correlated hopping Ti. 
Taking into account Ti leads to the concentration dependence of the hopping integrals and 
as a result to decreasing Uc- This effect shows itself the more the larger is the value of 
electron concentration n. A more detailed analysis of the correlated hopping Ti influence 
will be given in subsequent papers. 
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FIGURES 



FIG. 1. The possible configurations of the lattice sites. 

FIG. 2. Energy levels corresponding to the possible electron configurations of sites and the 
transitions between them. 

FIG. 3. The dependence of energy gap width AE/{U' — J) on the parameter ([/' — J)/{2w): 
the upper curve - {kT)/(2w) = 0.1; the middle curve - {kT)/{2w) = 0.05; the lower curve - 
{kT)/{2w) = 0. 

FIG. 4. The dependence of energy gap width AE/{U' — J) on the parameter {kT)/{2w): the 
upper curve - {U' - J)/{2w) = 0.74; the lower curve - {U' - J)/{2w) = 0.72. 

FIG. 5. The dependence of energy gap width AE/U on the parameter U/{2w) for n = 1 at 
different values of the parameter k = t[j/tij: the lower curve -k = 0; the middle curve -k = 0.2; 
the upper curve -k = 0.6. 

FIG. 6. The dependence of energy gap width AE/U on the parameter U/{2w) for n = 2 at 
different values of the parameter k = t^j/tif. the lower curve -k = 0; the middle curve -n = 0.1; 
the upper curve -k = 0.2. 

FIG. 7. The dependence of energy gap width AE/U on the parameter U/{2w) for n = 3 at 
different values of the parameter k = t'^j/Uj: the lower curve -n = 0; the middle curve -k = 0.1; 
the upper curve -k = 0.15. 

FIG. 8. The dependence of energy gap width AE/U on the parameter U/{2w) in the absence 
of correlated hopping (t[^ = *k = 0)- the lower curve ^ = 1, n = 3; the upper curve = 2. 

FIG. 9. The electron vs. interaction phase diagram showing the paramagnetic metal (PM) and 
paramagnetic insulator (PI) in absence of correlated hopping. 
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FIG. 10. The dependence of critical value {U/2w)c on the parameter of correlated hopping 
K = t'^Jtij-. the curve 1 - n = 1; the curve 2 - n = 2; the curve 3 - n = 3. 
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